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I.  INTRODUCTION 


Because  the  nosetip  shape  remains  symmetric  and  unchanged  throughout 
the  trajectory,  aerodynamic  uncertainties  for  reentry  vehicles  with  transpi¬ 
ration-cooled  nosetips  are  considerably  reduced.  Present  transpiration- 
cooled  nosetip  design  shapes  are  sphere -cones .  It  is  natural  to  ask  if 
additional  improvements  in  reentry  vehicle  performance  could  be  achieved 
through  the  selection  of  alternate  transpiration-cooled  nosetip  shapes.  One 
such  improvement  might  be  a  decrease  in  the  amount  of  coolant  required  by 
the  nosetip.  In  order  to  realize  such  a  decrease,  the-totaj.  integrated  heat 
transfer  to  the  nosetip  throughout  a  trajectory  must  be  reduced. 

The  objective  of  this  work  has  been  to  realistically  examine  the  poten¬ 
tial  for  decreasing  reentry  vehicle  nosetip  total  heat  transfer  through  the  use 

1-3 

of  alternate  nosetip  shapes.  In  earlier  work,  Baker  and  Kramer  derived 
analytic  expressions  for  nosetip  and  reentry  vehicle  total  trajectory  heat 
transfer.  In  Refs.  1  and  2  the  dependence  of  total  trajectory  heat  transfer  on 
reentry  parameters,  nosetip  or  body  scale,  and  geometry  (shape)  was  derived 
and  discussed.  Then,  in  Ref.  3,  this  information  was  utilized  and  combined 
with  classical  variational  calculus  procedures  and  Newtonian  pressure  and 
velocity  approximations  to  determine  minimum  heat  transfer  nosetip  shapes. 


Baker,  R.  L.  and  R.F.  Kramer,  "Evaluation  of  Total  Heat  Transfer  in 
Hypersonic  Flow  Environments,"  Report  No.  TR-0077{3550- 15)-3,  The 
Aerospace  Corporation,  El  Segundo,  California,  August  1978, 

^Baker,  R.  L.  and  R.  F.  Kramer,  "Reentry  Vehicle  Total  Heat  Transfer  for 
Fixed  Weight  and  Ballistic  Coefficient,  "  The  Aerospace  Corporation, 

El  Segundo,  California  (in  preparation). 

3 

Baker,  R.  L.  and  R.F.  Kramer,  "Nosetip  Shape  Optimization  for  Minimum 
Transpiration  Coolant  Requirements,  "  Progress  in  Aeronautics  and 
Astronautics .  5 9,  1978,  pp.  404-427. 


IRECEJXLliO  PMX  miffiD 


The  present  work  improves  upon  the  work  of  Ref.  3  by  using  inviscid 

•4  5 

flow  field  computer  codes  ’  to  compute  the  surface  pressure  distribution  in 
place  of  the  rough  Mewtonian  approximations.  The  heat  transfer  minimiza¬ 
tion  oroblem,  in  this  study,  is  treated  using  direct  optimization  methods  on 
flat-faced  geometries. 

In  Section  II,  the  required  total  heat  transfer  relationships  are  dis¬ 
cussed.  Then  in  Sections  III  and  IV,  the  optimization  problem  to  determine 
minimum  heat  transfer  nosetip  shapes  for  fixed  vehicle  ballistic  coefficient, 
nosetip  scale,  and  nosetip  fineness  ratio  is  defined,  and  the  present  results 
are  compared  with  earlier  approximate  results.  In  Section  V,  consideration 
is  given  to  the  use  of  sphere-cone,  minimum  heat  transfer,  and  truncated 
cone  nosetips,  with  fixed  fineness  ratio,  on  conical  reentry  vehicles  con¬ 
strained  to  have  specified  mass,  ballistic  coefficient,  and  volume.  These 
new  constraints,  which  take  into  account  the  interaction  of  nosetip  shape  with 
vehicle  performance  characteristics,  are  shown  to  favor  the  use  of  truncated 
cone  nosetips  for  achieving  even  lower  total  trajectory  nosetip  heat  transfer. 
No  attempt  was  made  to  optimize  nosetip  shape  to  achieve  a  minimum  under 
these  constraints. 


4 

Masson,  B.S.,  T.  D,  Taylor  and  R.  M.  Foster,  "Application  of  Godunov's 
Method  to  Blunt  Body  Calculations,  “  AIAA  Journal,  j7(4),  1969,  p.  6  94, 

^Baker,  R,  L.,  "Method  of  Characteristics  Computer  Program  Including 
Embedded  Shocks  and  Total  Enthalpy  Gradients  Normal  to  Streamlines,  " 
Report  No.  TOR-008 0 ( 55 50-0 1 ) -1 ,  The  Aerospace  Corporation,  San 
Bernardino,  California,  31  January  1980. 


-8- 


II*  HEAT  TRANSFER  PREDICTION  METHODS 


A.  INTEGRATED  HEAT  TRANSFER  EXPRESSIONS 

Wlien  expressions  for  the  laminar  ?.nd  turbulent  convective  heat- 

6  7 

transfer  coefficient  given  by  Vaglio-Eaurin  '  are  used  and  the  procedure 

g 

of  Allen  and  Eggers  is  followed,  the  rate  of  change  of  integrated  heat  trans¬ 
fer  with  altitude  dQ/dy  is  given  by 


dy 


-C 

Re 


Cp  :rr^^  r 

^co  CO  B  I 

J 

o  E, 


—  —  77  —  n  +  l 
/>  u  r  ds 
e  e  e 


■  Jo 


-I  n/n+i 


ds 


(1) 


where  n  =  1  for  laminar  flow  and  n  =  1/4  for  turbulent  flow.  When  an  inte¬ 
gration  in  nondimensiona  1  surface  distance  s  s  s/r  is  performed,  this 

max  ^  ’ 

expression  becomes 


o  E 


(2) 


where  the  integral  Iq 


is  given  by 


—  —  77  — n+1 

P  u  p  r  ds 
e  e  e 


(3) 


Vaglio- Laurin,  R.,  "Laminar  Heat  Transfer  on  Three-Dimensional  Blunt 
Nosed  Bodies  in  Hypersonic  Flow,  "  ARS  Journal,  29,  February  1959, 
pp.  123-129. 

^ Vaglio- Laurin,  R.,  "Turbulent  Heat  Transfer  on  Blunt  Nosed  Bodies  in  Two- 
Dimensional  and  General  Three-Dimensional  Hypersonic  Flow,  "  J,  Aero/ 
Space  Sciences ,  27,  January  19b0.  po.  27-36. 

8 

Allen,  H.  and  A.J,  Eggers,  Jr.,  "A  Study  of  the  Motion  and  A.erodynamic 
Heating  of  Ballistic  Missiles  Entering  the  Earth's  Atmosphere  at  High 
Supersonic  Speeds,  "  NACA  Report  1381.  1958. 


Allen  and  Eggers  integrated  an  expression  analogous  to  Eq.  (2)  representing 

and  as  functions  of  altitude  y  from  their  trajectory  analysis  and 

assuming  Re^  to  be  constant.  In  the  present  analysis,  the  shock- layer 

Reynolds  number  Re  =  [r  i/ZH  p  ]/u  has  also  been  represented  as  " 

’  o  max  ''  0*^0  o 

function  of  altitude  y. 

9 

Following  the  general  integration  procedure  of  Brunner  ,  the  total 
heat  transfer  for  an  entire  trajectory  Q  (Btu's)  is  obtained  by  integrating 
Eq.  (2),  The  result  for  laminar  flow  is  given  by 

K  ) 

and  for  turbulent  flow  by 

=  c'^f^  [/SlJ  ]  ^%[4/5.  (5] 


The  complete  derivation  of  Eqs,  (4)  and  (5)  is  given  in  Ref.  1.  The  ^  and 

2 

r  in  these  equations  are  the  vehicle  ballistic  coefficient  (lb  /ft  )  and 
max  m 

the  maximum  radial  coordinate  (ft)  of  the  nosetip,  respectively.  The  other 
parameters  are  given  by 


C 


L 


1.  866  x  10'^ 


[(y-H)/(y-i): 

(2w+3)A 


1/2 


C 


T 


4,  70 


X  10 

JTT 


-fa 


I  [(y+i)/(y-i)] 

(20/ +9) 


4/5 


(6a) 


(6b) 


9 

Brunner,  M.J.,  "Analysis  of  Aerodynamic  Heating  for  a  Re-entrant  Space 
Vehicle,"  J.  of  Heat  Transfer,  81,  August  1959,  pp.  223-229. 
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(7a) 


„l/2  ^3/2 
,1-  _  E  E 
*E  ‘ 


{sinBr 


1/59/5 

,T  _  ''e 

*E  "  ,  .  .  .1/5 


(sin^g) 


(7b) 


Q 


“  f  Pe“ 

Jo 


(8a) 


/: 


- - 5/4  — 

p  u  /i  r  d3 
*^£66 


(8b) 


B  -  [P°/(X;35in0j,)] 


=  [(2<i;+3)/4]B 


B^  =  [(2w+9)/10]B 


(9a) 

(9b) 

(9c) 


The  variables  in  these  equations  are  the  air  specific  heat  ratio  and  viscosity 
law  exponent  Y  and  u,  the  density  scale  height  factor  in  the  exponential 
atmosphere  X(ft'  ),  the  reentry  velocity  and  entry  flight  path  angle  Vj,(ft/sec) 

and  the  air  viscosity  evaluated  at  the  entry  stagnation  temperature 

^  o  3 

p_  (lb  /ft-sec),  the  ballistic  factor  B,  the  sea  level  air  density  p  (lb  /ft  ), 

jid  y  rrx 

and  the  nondimensional  heat  transfer  integrals  Iq  and  Iq  which  depend  upon  an 

integration  in  normalized  surface  running  length  of  the  local  nondimensional 

boundary  layer  edge  density  p^  =  velocity  u^  =  u^/v/2H^,  viscosity 

U  -  U  I U  ,  a-nd  local  normalized  radial  coordinate  t  =  r/r  .  The  functions 
^e  max 
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appearing  in  Eqs.  (4)  and  (5)  are  the  error  function  and  the  incomplete  gamma 

function,  respectively.  These  functions  are  given  in  terms  of  the  appropriate 

Li  T 

laminar  and  turbulent  modified  ballistic  factors  B  and  B  in  Fig.  1. 

It  is  now  evident  that  if  the  reentry  parameters  and  0^,  vehicle 

ballistic  coefficient  3,  and  maximum  nosetip  radial  coordinate  r  are 

fixed,  then  the  total  laminar  heat  transfer  Q  will  be  minimized  if  we  mini- 

L  T 

mize  the  integral  Iq.  Similarly,  in  the  turbulent  case,  Q  is  minimized  by 

minimizing  the  integral  Iq. 

B .  EVA  LUATION  OF  Iq  AND  Iq 

L  T 

The  heat  transfer  integrals  1^  and  Iq,  given  by  Eqs.  (8a)  and  (8b), 
respectively,  depend  only  upon  the  nosetip  geometry,  i.  e.,  upon  the  shape 
contour  of  the  nosetip.  The  basic  geometry  considered  in  this  and  in  earlier 
work  is  shown  in  the  inset  of  Fig.  2.  It  consists  of  a  nosetip  having  a  flat  face 
of  height  r(0),  to  be  determined,  a  maximum  radial  coordinate  at  the 
back  ^  smoothly  varying  contour  r  *  r(z)  from  the  flat  face  expan¬ 

sion  corner  to  the  back  of  the  nosetip.  The  flat  face  height  r(0)  is  not 
constrained  to  be  nonzero.  Thus,  the  minimum  heat  transfer  shape. deter¬ 
mined  is  not  required  to  have  a  flat  face.  Note  that  the  local  nosetip  coordi¬ 
nates  r  and  z  and  the  surface  running  length  s  are  related  to  the  local  nosetip 
surface  angle  Ct  by  the  following  simple  trigonometric  relationships: 

57  '  ■  if  = 

In  Ref.  3,  the  local  nondimensional  boundary  layer  edge  properties  p^, 

u^,  and  were  represented  as  'lunctions  of  Ct,  for  ct  <  90“,  by  the  simple 

hypersonic  Newtonian  approximations  for  pressure  and  velocity.  Over  the 

flat  face  portion  of  the  nosetip,  Ct  =  90°,  the  product  p^xi^fl^  was  assumed  to 

vary  linearly  with  the  normalized  distance  to  the  sonic  point  at  the  expansion 

corner.  In  the  present  work,  (.he  flow  pressure  distributions  were  obtained 

4  5 

numerically  from  accurate  inviscid  flow  field  calculation  procedures.  ’ 
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Fig.  1.  Laminar  and  Turbulent  Heat  Transfer  Functi 


A  schematic  diagram  illustrating  the  calculation  procedure  is  shown  in 

4 

Fig.  2.  The  time-dependent  (Godunov  )  computer  code  required  about  30  min 
of  CDC  7600  computer  run  time  to  calculate  the  blunt  body  portion  of  the  flat 
face  flow  field.  However,  it  was  only  necessary  to  perform  this  calculation 
once,  since  this  portion  of  the  flow  field  does  not  depend  upon 
upon  the  shape  contour  after  the  expansion  corner.  The  supersonic  portion 
of  the  flow  field  in  all  cases  was  calculated  using  a  fast  running  rotational 
characteristics  computer  program^  which  required  only  two  to  three  sec  of 
computer  time  per  specified  flat  face  nosetip  shape. 

The  normalized  product  can  be  written  in  terms  of  the  normalized 

pressure  p^  =  Pg/p^  as 


P 

p  fj  =  _£_£ 

^e  e  p  fi 
*0  0 


_  [(y-i)tj+i]/y 

Pe 


(11) 


Also,  assuming  an  isentropic  expansion,  the  local  nondimensional  velocity  u^ 
can  be  written  in  terms  of  p  as 


Thus,  Iq  is  a  function  of  the  specific  heat  ratio  y,  the  temperature  exponent 

in  the  viscosity  law  <*>,  and  body  geometry.  In  this  work,  the  hypersonic 

values  y  =  1,2  and  <J  =  1.0  have  been  used  throughout.  The  integrals  li*  and 
'p  w 

Iq  were  calculated  numerically  from  Eqs.  (8a)  and  (8b)  for  specified  geometry 

using  the  calculated  surface  pressure  distribution  for  that  geometry  and 

Eqs.  (11)  and  (12). 
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Flowfield  Calculation  Methods 


m.  PROBl^M  DEFINITION 


A.  MATHEMATICAL  FORMULATION 

In  this  section,  we  will  be  concerned  with  the  determination  of  those 
flat-faced  shapes  that  produce  minimum  values  of  Iq.  Nosetip  shapes  with 
a  flat  forward  face  of  unspecified  height  r  =  r(0)  will  be  considered  (see 
Fig.  2).  Normalizing  all  coordinates  by  define 


=  z/] 


max 


y  - 


max 


(13) 


Following  the  general  procedures  used  in  Ref.  3,  the  Iq  integrals  given 
by  Eqs.  (8a)  and  (8b)  may  be  expressed  as 


*  ^^Q^face 


nose 


{14a) 


=  n^I  + 


/ 


2/7r  —  -  -  n+1 

p  M  H  y 
^e  e  e 

cos  d 


dx 


(X4j) 


where  n  a  1  for  laminar  flow  and  n  =  1/4  for  turbulent  flow.  The  fineness 

ratio  T=  2r  / L  is  to  be  specified, 
max  ^ 

The  value  of  the  parameter  G  defines  the  heat  transfer  to  the  flat  face 
portion  of  the  nosetip.  The  G  for  both  laminar  and  turbulent  flow  was  deter¬ 
mined  in  the  present  work  using  calculated  p  values  from  the  finite  element 

4  ®  _  _  _ 

(time -dependent)  flow  field  computer  code  ,  evaluating  and  u^  from  Eqs. 

(11)  and  (12),  respectively,  and  numerically  integrating  Eqs,  (8a)  and  {8b). 
The  calculated  values  of  G  are 


=  0.0375 


G  =  0.  0449 


(15a) 

(15b) 
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Substituting  Eqs.  (11)  and  (12)  into  Eq.  {14b),  the  fiinctional  to  be 
minimized  may  be  written  in  the  following  form; 


rZ/r 

Iq  *  Pg.  y)dx  (16) 

Jo 


Because  of  the  coupling  with  the  method-of-characteristics  computer  program 
for  evaluation  of  the  integral  in  Eq.  (16),  it  was  decided  not  to  attempt  to 
derive  an  Euler  equation  corresponding  to  this  variational  problem  as  done 
in  Ref.  3.  Instead,  a  purely  numerical  minimization  approach  was  used. 

The  shape  y(x)  was  approximated  by  a  truncated  Fourier  series  which  satis¬ 
fies  the  boundary  conditions  y(0)  =  dy/dx(0)  =  b  (constant  b  specified  in  Sec¬ 
tion  III-B)  and  y(2/T)  =  1. 


y(x) 


N 


+  bx  + 


Z“i 

i=l 


cos 


2i-l 


ITTK 

2 


(17) 


A  nosetip  shape  is  thus  defined  by  specifying  the  N  parameters 
flL  (i«l,  .  .  .  ,  N).  For  a  given  shape,  the  flow  field  may  be  found  using  the 
method  of  characteristics  code.  This  determines  p^,  u^,  and  along  the 
body,  through  Eqs.  (11)  and  (12),  and  the  integral  in  Eq.  (16)  may  then  be 
evaluated.  The  variational  problem  has  thus  been  replaced  by  the  finite  para¬ 
meter  minimization  problem;  minimize 

Iq  (ttj - ,aj^).  (18) 

This  problem  was  solved  iteratively  using  the  Davidon-Fletcher-Powell 
Method.  ^  It  was  determined  that  N  =  3  produced  sufficient  accuracy  for  the 
present  calculations. 

^^Greenstadt,  J.,  “Variations  on  Variable- Metric  Methods,  "  Math,  of  Comp. 
1970,  pp.  1-22.  - 
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B. 


TRANSVERSAUTY  CONDITION 


Since  y(0)  is  not  specified,  the  minimization  of  the  functional  Iq,  as 
given  by  Eq.  (16),  is  a  variable  end  point  problem  of  the  calculus  of  varia¬ 
tions,  and  the  following  trans versality  condition  must  be  satisfied  at  x  =  0. 

^  -  G  y{0)’'^^  =  0  (19) 

dy 


The  flow  downstream  from  the  expansion  corner  is  supersonic,  there¬ 
fore,  conditions  at  x  =  0  cannot  be  influenced  by  conditions  downstream. 

Hence,  the  transversality  condition  Eq.  (19)  may  be  evaluated  using  the 

12 

Prandtl-Meyer  relations  to  represent  the  flow  variables  as  a  function  of 
flow  angle  through  the  corner  expansion.  This  leads  to  the  following  expres¬ 
sion  for  dT /Qy, 


dF 

dy 


^{y-i)6i+i]/y 

^l^e 


+  K- 


L(y-l)a;+l]/y 
^e  • 


(y-D/y 


where 


(20) 


Kj  =  sinOI  + 


(y-i)w+i 

yM^ 

y 

.  '/m^-i  . 

cosa 


(21a) 


K,  =  -sina  + 


1 

yM^ 

y 

cosa 


(21b) 


and  M  is  the  Mach  number  corresponding  to  the  flow  angle  a  in  the  Prandtl- 
Meyer  expansion. 

Miele,  A.  (ed.  ),  '•Theory  of  Optimum  Aerodynamic  Shapes,  '•  Applied 
Mathematics  and  Mechanics.  Academic  Press,  N.  Y.,  1965. 

Liepmann,  H.  W.  and  A.  Roshko,  Elements  of  Gasdynamies,  John  Wiley  & 
Sons,  Inc.,  N.Y.  (1957),  pp.  92-100. 
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If  y(0)  *  0,  Eq.  (19)  may  be  divided  by  y{0)^^^.  Thereafter,  substi¬ 
tuting  for  5F/9y  fi'om  Eq.  (2U),  the  resulting  equation  is  a  nonlinear  equation 
in  y(0).  The  numerical  solution  of  this  equation  yields  the  following  values: 

y(0)^  =  0.  3716  (22a) 

y(0)^  =  0,  314/  (22b) 

Equations  (22a)  and  (22b)  fix  the  slope  of  the  minimum  heat  transfer 
nosetip  shapes  at  x  =  0  for  the  laminar  and  turbulent  cases,  respectively. 

The  angles  corresponding  to  these  slopes  are  20,39°  and  17.47°,  respectively. 
The  values  of  y(0)  at  x  =  0,  i.  e.,  the  flat  face  heights,  are  determined  itera¬ 
tively  as  part  of  the  overall  solution  procedure  described  in  the  problem 
formulation  discussion  above. 
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IV.  RESULTS  AND  DISCUSSION 


Results  obtained  by  use  of  the  direct  optimization  methods  discussed 
above  are  presented  herein  for  the  cases  of  both  laminar  and  turbulent  flow  in 
the  boundary  layer.  In  all  cases,  the  specific  heat  ratio  y  and  the  ^hscosity 
law  exponent  <*>  have  been  assumed  to  be  the  hypersonic  values  1,  2  and  1.  0, 
respectively. 

A.  MINIMUM  HEAT  TRANSFER  NOSETIP  SHAPES 

Calculated  minimum  laminar  and  turbulent  heat  transfer  shapes  for 
fineness  ratios,  T  =  2r  / L,  of  1,2,4  and  8  are  shown  in  Fig.  3.  The  flat 
face  height  for  the  minimum  heat  transfer  shape,  determined  in  each  case  as 
part  of  the  solution,  decreases  quite  rapidly  from  a  large  fraction  of 
for  T  =  8  to  35%  of  r  for  T  =  1.  In  all  cases,  the  r(0)  is  greater  for  the 
turbulent  flow  solution  than  for  the  laminar  case  with  the  same  T.  This  is 
probably  due  to  the  smaller  exponent  of  r(0)  in  the  expression  for  [iQlfagg- 
Eq.  (14b),  in  the  case  of  turbulent  flow  than  in  the  case  of  laminar  flow. 

For  comparison,  the  minimum  heat  transfer  shapes  calculated  earlier 
in  Ref.  3  using  the  Newtonian  approximations  are  also  shown  in  Fig.  3. 

In  all  cases,  the  flat  face  heights  are  greater  than  those  for  the  present  solu¬ 
tions.  This  is  especially  true  for  the  T  =  1  case.  The  main  reason  for  this  is 
the  larger  flow  angles  at  the  expansion  corner  (17-20°  instead  of  3-4°), 
indicated  by  the  present  trans versality  condition  using  the  Prandtl-Meyer 
equations  in  place  of  the  Newtonian  approximations.  In  the  Newtonian  ap¬ 
proximation  solutions,  so  little  difference  is  predicted  between  the  minimum 
heat  transfer  laminar  and  turbulent  shapes  for  these  T  values  that  only  one 
shape  is  shown.  The  Newtonian  approximation  solutions  each  have  monotoni- 
catly  decreasing  slopes  after  the  expansion  corner,  whereas  the  present 
minimum  heat  transfer  nosetip  shapes  are  inflected  for  t2;  2.  A  truncated 
sphere  or  a  5°  sphere-cone  shape  is  shown  in  each  case  for  shape  comparison 
purposes. 
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Minimum  Heat  Transfer  Nosetip  Shapes 


B. 


REL.ATIVE  HEAT  TRANSFER  RATES 


Calculated  minimum  laminar  and  turbulent  values  for  the  shapes 
shown  in  Fig.  3  are  given  in  Fig.  4.  The  earlier  results,  also  shown,  had 
indicated  that  the  laminar  and  turbulent  Iq  values  for  a  minimum  heat  trans¬ 
fer  nosetip  shape  were  up  to  three  to  four  times  smaller  than  those  for  a 
truncated  sphere  having  the  same  fineness  ratio.  The  present  results  indicate 
calculated  Iq  values  for  truncated  spheres  much  lower  than  those  of  the 
earlier  results  and  calculated  minimum  Iq  values  higher  than  those  indicated 
from  the  approximate  calculations.  Overall,  this  dictates  that  the  potential 
reduction  in  heat  transfer,  obtained  by  employing  a  minimum  heat  transfer 
nosetip  shape  instead  of  a  sphere-cone  shape,  is  m.uch  less  than  predicted 
earlier. 

The  reasons  for  this  behavior  are  related  to  the  type  of  errors  inherent 

in  the  Newtonian  approximations.  In  general,  the  Newtonian  velocity 

u  =  cosOt  is  as  much  as  a  factor  of  2  too  high,  whereas  the  Newtonian 
®  —  2 

pressure  p^  =  sin  <X  is  reasonably  accurate  on  the  blunt  portion  of  the  nose¬ 
tip  (C*  >  45)  and  becomes  even  more  than  a  factor  of  2  too  low  for  very  small 
body  angles  (d  <  5).  In  the  results  shown  in  Fig.  4,  the  velocity  error 
tends  to  dominate  in  the  spherical  shape  calculations,  whereas  the  pressure 
error  dominates  in  the  minimum  heat  transfer  shape  predictions. 

A  comparison  of  the  calculated  numerical  values  of  Iq  for  tie  present 
minimum  heat  transfer  shapes  to  those  for  truncated  spheres  or  sphere- 
cones  having  the  same  fineness  ratio  T  is  shown  in  Table  1.  The  maximum 
reduction  in  Iq  is  about  35%.  As  a  m.atter  of  interest  in  later  discussion, 
the  values  of  the  drag  integral  Ij^  for  the  various  shapes  are  also  shown. 
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Fig,  4.  Laniinar  an<i  Turbulent  Values 


Table  i.  Calculated  Hea 
Drag  Integrals 


V.  REENTRY  VEHICLE  APPLICATIONS 


The  minimum  heat  transfer  shapes  determined  in  the  previous  sections 
were  obtained  assuming  that  the  reentry  conditions,  ballistic  coefficient 
and  nosetip  scale  r  could  be  fixed  independently.  If  we  consider  the 
interaction  of  yS  and  then  these  quantities  can  no  longer  be  specified 

independently.  In  this  section,  sphere-cone  reentry  vehicles  constrained 
to  have  constant  mass  m  and  ballistic  coefficient  p  and  specified  cone  angle 
and  nosetip  fineness  ratio  r  are  first  considered.  Then,  fixed  cone  half-angle 
reentry  vehicles  having  minimum  heat  transfer  or  truncated  cone  nosetips 
with  specified  T  are  considered  and  constrained  to  have  constant  m,  and 
volume  V.  In  all  cases  considered,  the  nosetip  and  reentry  vehicle  geometry 
is  specified  a  priori,  the  constraints  are  applied,  but  no  optimization  pro¬ 
cedure  is  attempted.  This  determines  the  values  or  r„  and  r  and.  sub- 

B  max  ' 

sequently,  Q. 

L  3/2 

Since  by  Eq.  (4)  Q  is  proportional  to  and  according  to  Eq.  (5) 

Q  is  proportional  to  r^^  ,  the  total  heat  transfer  Q  may  also  be  reduced 

max  ' 

by  decreasing  r  .  The  reduction  in  this  case  is  due  to  the  decrease  in 
surface  area  of  the  nosetip. 

Reduction  in  Q  by  reducing  is  illustrated  for  the  selected  nosetip 

shapes  with  the  indicated  constraints  in  the  following  subsections. 

A.  SPHERE-CONE  NOSETIPS.  CONSTANT  m  AND  0 

To  illustrate  the  reduction  in  Q  by  means  of  decreasing  r  ,  consider 

®  max 

the  following  baseline  sphere -cone  reentry  vehicle: 

m  s  400  lb 

m 

B  =  2000  lb  /ft^ 
m 

0=7° 

c 
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IMBCSnUO  Fica  BLAMUNOV  flUflO 


r.,  =  Z.  56  in.,  r  *2.61 
N  max 


T  =  2  r  /L  =  2.  0. 

nosetip  max 


Our  present  objective  is  to  reduce  r  so  as  to  reduce  Q,  but  to  do  so  main- 

max  g  22 

taining  m  and  constant.  Numerically  calculated  results  '  satisfying  the 
required  constraints  are  shown  in  Fig.  5  and  its  inserted  table.  The  reentry 
vehicle  base  radius  and  mass  r^  and  m,  the  ballistic  coefficient  and  a 
drag  integral  are  related  by 


m 


2  , 


{23a) 


(23b) 


where 


■/ 


p  sinOfr  ds 


(24a) 


{24b) 


Reducing  while  maintaining  9^  constant  decreases  the  drag 

integral  Requiring  m  and  0  to  be  constant  thus  means,  by  Eq.  (23b),  that 

the  vehicle  base  radius  must  increase.  This  is  the  behavior  seen  for  the 
three  reentry  vehicles  shown  in  Fig,  5  and  the  corresponding  vehicle  para¬ 
meters  given  in  the  inserted  table. 


13 


Inouye,  M,,  J.  V.  Rakich  and  H.  Lomax,  "A  Description  of  Numerical 
Methods  and  Computer  Programs  for  Two-Dimensional  and  Axisymmetric 
Supersonic  Flow  over  Blunt-Nosed  and  Flared  Bodies,  ”  NASA  TN-D-2970, 
August  1965. 
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VEHICLE 


Since  all  other  factors  are  constant,  the  nosetip  total  heat  transfer  Q 
for  these  three  vehicles  is  simply  proportional  to  a  power  of  i.e. 


/l. 83^/2 

[zruj 


0.687; 


0,  528 


0.849; 


0.  822 


The  maximum  reduction  in  nosetip  total  heat  transfer  of  47.2%  is  obtained  at 
the  expense  of  an  enormous  increase  in  the  volume  and  surface  area  of  the 
reentry  vehicle.  In  many  instances,  such  increases  will  not  be  tolerable  in 
view  of  volume  and  heat  shield  weight  requirements, 

B.  MINIMUM  HEAT  TRANSFER  NOSETIPS. 

CONSTANT  m.  B.  AND  V 

The  interaction  of  nosetip  shape  with  reentry  vehicle  characteristics 
and  performance  is  considered  below  for  the  case  of  minimum  heat  transfer 
nosetips  on  conical  reentry  vehicles.  Since  volume  increases,  such  as 
those  encountered  in  the  sphere-cone  geometry  results  just  discussed,  are 
undesirable,  this  additional  constraint  will  also  now  be  considered. 

The  volume  of  an  axisymmetric  reentry  vehicle  can  be  defined  in 

2 

terms  of  a  volvime  integral  I,^^  and  the  base  radius  r^  by 
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where  ly  is  given  by 


and  from  this  we  can  write 


(25b) 


(25c) 


If  we  require  a  reentry  vehicle  to  simultaneously  satisfy  a  ballistic  coefficient 
constraint  and  a  volume  constraint,  then  equating  the  base  radius  expressions 
given  by  Eqs.  (23b)  and  (25c)  gives  the  following  required  relationship  between 
m.  V.  and  ly 


0  = 


=  2.42 


ll/2 

D 

-pr 

r 


(26) 


The  calculated  value  of  this  dimensionless  parameter  for  the  baseline 
vehicle  considered  in  Subsection  V-A  is  0.  293.  Our  present  objective  is  to 
determine  the  nosetip  total  heat  transfer  for  a  T  =  2  minimum  turbulent  heat 
transfer  nosetip  on  a  7°  conical  reentry  vehicle  for  which  the  m,  and  V 
are  the  same  as  for  the  baseline  sphere-cone  vehicle. 

For  specified  nosetip  and  conical  afterbody  geometry,  0  can  easily  be 

Li  T 

determined  in  the  same  manner  as  the  integrals  Iq  andT^,  i.e.,  using 
numerically  calculated  values  of  p^  to  evaluate  1^^  as  given  by  Eq.  (24a)  and 
Eq.  (25n)  to  evaluate  ly.  When  this  is  done  for  the  present  geometry,  it  is 
found  that  for 

r 

0  =  0.293;  — — -  =  0.260  (27) 

The  corresponding  value  of  r  /r_  for  the  baseline  sphere-cone  vehicle  is 
0. 287. 
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Designating  the  sphere -cone  reentry  vehicle  as  1  and  the  minimum 
heat  transfer  nosetip  vehicle  as  2,  we  also  have  for  fixed  and  m  from 
Eq.  {23b)  and  the  calculated  integrals 


1/2 

1 

TTT 


0.998 


(28) 


Thus,  combining  these  relationships 


’max*  2  ^^max'^^B*^ 


dr 


(r 


max 


^"‘max^'^B*^  (Ij^) 


1 

TTT 


=  0. 904 


(2  9) 


and,  finally,  making  use  of  Eq.  (5)  and  Table  1 


Q 

Q 


’  (r  )  " 

max  2 

9/5 

T 

2 

(r  ) 

L  max  j  . 

,rTv 

4/5 


0.698 


(30) 


Thus,  the  nosetip  total  turbulent  heat  transfer  is  reduced  30.2%,  even  though 
the  reentry  vehicle  volume  is  held  constant  along  with  the  mass  and  ballistic 
coefficient.  The  main  reason  this  can  be  done  is  the  intrinsically  higher  drag 
coefficients  for  reentry  vehicles  with  minimum  heat  transfer  nosetips. 

The  basic  quantity,  in  addition  to  I_,  needed  to  compare  nosetip  (T  =  2) 

■  1  /? 

total  heat  transfer  as  inEqs.  (29)  and  (30)  is  (r  /r^]/{l^)  '  .  This  para- 

0  0  o  mftx 

meter  is  shown  for  5,7,  and  9  conical  reentry  vehicles  with  T=  2  minimum 
turbulent  heat  transfer  nosetips  in  Fig.  6.  The  corresponding  results  for  5°, 

7  ,  and  9°  conical  reentry  vehicles  with  T  =  2  sphere-cone  nosetips  are  also 
shown.  The  baseline  vehicle  and  the  minimum  heat  transfer  nosetip  vehicle 
discussed  above  are  indicated  by  an  "x"  on  the  7^  curves. 
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Using  Fig.  6,  comparisons  such  as  that  made  above  can  easily  be  made 
for  a  wide  range  of  reentry  vehicle  parameters.  The  correlating  parameter 
is  0  a 

C.  TRUNCATED-CONE  NOSETIPS. 

CONSTANT  m,  B  ,  AND  V 

From  the  above  results,  the  desirability  of  a  "high  drag"  nosetip  for  the 
present  puTpose  of  reducing  nosetip  total  heat  transfer  can  be  seen.  An 
obvious  new  objective  thus  becomes  finding  a  nosetip  shape  which  has  large 
drag  (high  Ij^)  and  low  heat  transfer  (small  Iq).  Comparison  of  these  quantities 
for  truncated  cone  nosetips,  minimum  heat  transfer  nosetips,  and  sphere- 
cone  nosetips  is  shown  in  Table  2.  Comparison  of  corresponding  quantities 
for  minimum  heat  transfer  nosetip  shapes  and  truncated  cones  indicates  that 
the  latter  Iq  values  are  only  2  to  3%  greater,  whereas  the  latter  values 
are  28-47%  greater  for  T>  2. 

This  means  that  using  a  truncated  cone  nosetip  in  place  of  a  sphere- 
cone  nosetip  and  maintaining  constant  m,  ;S,  and  V,  as  before,  could  potentially 
give  even  larger  reductions  in  Q  than  using  the  previously  determined  minimum 
heat  transfer  nosetip  shapes.  The  quantity  (r  /r_)/I^^^  for  5°,  7°,  and 
9  conical  reentry  vehicles  with  5  ,  t  ,  and  9  truncated  cone  nosetips, 
respectively,  is  shown  as  a  function  of  0  in  Fig.  7. 

For  the  7°  vehicles  marked  by  the  small  "x"  on  the  7°  curves,  the 
truncated  cone  nosetip  total  turbulent  heat  transfer  is  57.  9%  of  that  for  the 
reference  7°  sphere -cone  vehicle  having  the  same  m,  yS,  and  V.  Recall  that 
the  corresponding  reduction  employing  a  minimum  heat  transfer  nosetip  was 
to  69.  8%  of  the  sphere-cone  nosetip  value.  These  results  do  not  negate  the 
minimum  heat  transfer  shape  results  determined  earlier.  The  minimum 
heat  transfer  shapes  minimize  0(1,-.)  for  fixed  B  and  r  as  seen  in  Table  2. 
For  fixed  m,  ,  and  V,  truncated  cone  nosetips  allow  a  greater  reduction  in 
^max  they  have  larger  values);  this  results  in  a  greater  reduction 

in  Q  by  Eq.  (5)  even  though  Iq  for  these  nosetips  is  slightly  larger  than  for 
the  minimum  turbulent  heat  transfer  shapes. 
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Calculated  Nosetip  Heat  Transfer  and 
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Fig.  7.  Nosetip  Scaling  Paranneter  for  Sphere-cone  and 
Truncated  Cone  T=  Z  Nosetios 


Many  conaparisons  can  be  made  using  the  results  given  in  Figs,  6  and  7 
in  the  manner  illustrated  herein.  In  general,  the  magnitude  of  the  potential 
reduction  in  Q  appears  to  achieve  a  maximum  at  a  value  of  ^  which  increases 
as  the  cone  half-angle  increases.  The  maximum  reduction  for  truncated 
cone  nosetips  (T  =  2)  approaches  50%.  It  is  very  likely  that  reductions  greater 
than  this  could  be  achieved  by  simultaneously  optimizing  both  the  nosetip  and 
the  reentry  vehicle  shape  with  full  consideration  of  their  interaction. 

Fixed-shape  flat-faced  nosetips,  such  as  those  considered  in  this  study, 

may  possess  other  advantages  for  reentry  vehicle  applications  in  addition  to 

reducing  nosetip  total  heat  transfer.  They  may  significantly  reduce  the  nose- 

1  4 

tip  boundary  layer  transition  altitude.  In  addition,  the  rapid  recompression 
occurring  on  the  nosetip  surface,  just  after  the  overexpansion  at  the  flat- 
face  corner,  could  possibly  serve  as  a  very  effective  adverse  pressure  gradient 
boundary  layer  trip  for  frustum  transition.  This  might  assure  uniform  move¬ 
ment  of  the  frustum  transition  front  location  to  the  nosetip  region.  This 
becomes  potentially  ev^en  more  attractive  when  coupled  with  significant  re¬ 
duction  in  the  boundary  layer  unit  Reynolds  nunnber  on  the  rear  portion  of 
the  reentry  vehicle  frustum  due  to  the  flat- nose  shape. 


Auerbach,  I.,  "Influence  of  Nosetip  Shape  on  Boundary  Layer  Transition  in 
Arc  Heated  Jets,  "  AIAA  Paper  78-235,  AIAA  i6th  Aerospace  Sciences 
Meeting,  Huntsville,  Ala.,  January  1978. 
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VI.  SUMiVLARY  AND  CONCLUSIONS 


Flat-face  nosetip  shape  contours  that  minimize  reentry  trajectory 
total  nosetip  heat  trans fe r  for  fixed  vehicle  ballistic  coefficient,  nosetip 
fineness  ratio,  and  scale  have  been  determined  using  a  direct  optimization 
approach.  The  method  utiEzed  relies  upon  analytic  expressions  for  the 
total  integrated  heat  transfer  to  a  nosetip  for  a  constant  ballistic  coefficient 
reentry  trajectory.  These  expressions  are  based  upon  the  use  of  boundary 
layer  integral  heat  transfer  prediction  methods.  The  surface  pressure  distri¬ 
butions,  required  to  evaluate  the  heat  transfer  and  drag  functions,  have  been 
determined  numerically  from  accurate  inviscid  flow  field  computer  codes. 

Minimum  heat  transfer  flat-faced  shapes  have  been  found  for  fineness 
ratios  of  1,  2,  4,  and  8.  The  flat- face  height,  for  both  the  laminar  and  tur¬ 
bulent  solutions,  decreases  as  the  fineness  ratio  decreases.  These  shapes 
are  qualitatively  similar  to  those  obtained  earlier  using  Newtonian  flow  field 
approximations.  However,  the  potential  reduction  in  total  heat  transfer 
through  the  use  of  minimum  heat  transfer  nosetip  shapes  in  place  of  the 
normal  sphere-cone  geometry  is  much  less  than  predicted  earlier.  The 
maxinrium  heat  transfer  reduction  due  to  nosetip  shaping  alone,  i.e.,  fixed 
nosetip  scale,  is  about  30%. 

The  optimization  problem  solved  to  determine  the  minimum  heat  trans¬ 
fer  shapes  neglects  the  interaction  of  the  nosetip  shape  with  reentry  vehicle 
performance  characteristics.  The  use  of  minimum  heat  transfer  and 
truncated  cone  nosetips  on  conical  reentry  vehicles  with  this  interaction 
considered  has  also  been  investigated.  Because  these  shapes  are  high  drag- 
low  heat  transfer  configurations,  additional  decreases  in  nosetip  total  heat 
transfer  can  be  realized  due  to  the  reduced  scale  of  the  required  nosetip  in 
comparison  to  that  for  a  sphere-cone  nosetip  on  a  reentry  vehicle  having  the 
same  mass,  ballistic  coefficient,  and  volume.  No  attempt  to  optimize  the 
nosetip  shape  for  these  additional  constraints  was  made. 
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Total  heat  transfer  reductions,  with  these  new  constraints,  of  40%  are 
possible  using  the  previously  defined  T=  2  minimum  heat  transfer  nosetip 
shapes,  and  reductions  approaching  50%  are  possible  using  the  higher  drag 
T  s  2  truncated  cone  nosetip  shapes.  Scaling  relationships  have  been  derived 
and  presented  v/hich  allow  the  total  trajectory  nosetip  heat  transfer  for  these 
nosetip  shapes  to  be  easily  related  to  that  for  sphere-cone  nosetips  on  conical 
reentry  vehicles  constrained  to  have  the  same  mass,  ballistic  coefficient, 
and  volume. 


